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ABSTRACT 

With frequent flaring activity of its relativistic jets, Cygnus X-3 (Cyg X-3) is one 
of the most active microquasars and is the only Galactic black hole candidate with 
confirmed high energy y-ray emission, thanks to detections by Fermi / LAT and AGILE. 
In 2011, Cyg X-3 was observed to transit to a soft X-ray state, which is known to 
be associated with high-energy y-ray emission. We present the results of a multi- 
wavelength campaign covering a quenched state, when radio emission from Cyg X-3 
is at its weakest and the X-ray spectrum is very soft. A giant (~ 20 Jy) optically thin 
radio flare marks the end of the quenched state, accompanied by rising non-thermal 
hard X-rays. Fermi / LAT observations (E ^100 MeV) reveal renewed y-ray activity 
associated with this giant radio flare, suggesting a common origin for all non-thermal 
components. In addition, current observations unambiguously show that the y-ray 
emission is not exclusively related to the rare giant radio flares. A 3-week period of 
y-ray emission is also detected when Cyg X-3 was weakly flaring in radio, right before 
transition to the radio quenched state. No y rays are observed during the ~ one-month 
long quenched state, when the radio flux is weakest. Our results suggest transitions 
into and out of the ultrasoft X-ray (radio quenched) state trigger y-ray emission, 
implying a connection to the accretion process, and also that the y-ray activity is 
related to the level of radio flux (and possibly shock formation), strengthening the 
connection to the relativistic jets. 

Key words: black hole physics — X-rays: binaries - gamma-rays: star — ISM: jets 
and outflows — radio continuum: stars — stars: individual (Cyg X-3) 
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1 INTRODUCTION 

Galactic and extra-galactic accreting systems containing a 

neutron star or a black hole can produce outflows containing 
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energetic particles that are accelerated away from the com- 
pact object up to relativistic speeds in collimated jets. These 
high energy particles, entangled in the jet magnetic field, 
lose their energy via synchrotron, inverse Compton emission 
and/or adiabatic losses, or via pion production in the case of 
baryonic jets, resulting in a broad-band spectrum from radio 
up to high energy 7 rays. (e.g. Atoyan & Aharonian 1999; 
Georganopoulos, Aharonian & Kirk 2002; Romero et al. 
2003). 

Cyg X-3 was one of the first sources to be discov- 
ered in the early days of X-ray astronomy (Giacconi et al. 
1967). It is a system consisting of a Wolf-Rayet star 
(van Kerkwijk et al. 1992) and a compact object (most 
likely a black hole). With a short 4.8-hr orbital period 
(Parsignault et al. 1972), Cyg X-3 lies at a distance of the 
order of ~ 7 kpc (Ling, Zhang & Tang 2009). The X-ray 
spectrum from Cyg X-3 changes between hard and soft 
states akin to those observed in other accreting X-ray bi- 
naries, and is heavily absorbed at low energies by the inter- 
vening dense stellar wind (Szostek, Zdziarski & McCollough 
2008; Hjalmarsdotter et al. 2009). Cyg X-3 is also known 
for the recurrent activity of its relativistic jets that make 
it one of the brightest Galactic transient radio sources (e.g. 
Mioduszewski et al. 2001; Miller- Jones et al. 2004). 

AGILE and Fermi / LAT reported concurrent detections 
of Cyg X-3 in high energy 7 rays (>100 MeV), closely related 
to the activity of the relativistic jet during the soft X-ray 
state (Tavani et al. 2009; Fermi LAT Collaboration et al. 
2009, hereafter FLC09). The 7 -ray emission measured by 
the Fermi / LAT was found to be modulated on the orbital 
period, securing the identification (FLC09) . A very short 7 - 
ray flare was reported later during a short transient softening 
of Cyg X-3 in 2010 (Corbel & Hays 2010; Bulgarelli et al. 
2010; Williams et al. 2011). No evidence for emission above 
250 GeV has been found by MAGIC (Aleksic et al. 2010). 

Cyg X-3 is the first binary hosting an accreting com- 
pact object and relativistic jet (a.k.a. microquasar) to be 
detected in 7 rays. It is clearly an accreting source unlike 
gamma-ray binaries that are more naturally explained by 
pulsar spindown power (Dubus 2006). The nature and lo- 
cation of the non-thermal processes that bring particles to 
high energies, the relationship between these processes, jet 
launching and accretion state are long-standing questions for 
both microquasars and active galactic nuclei that stand to 
benefit from this detection. In Cyg X-3, jet emission can be 
followed and resolved spatially in radio, the accretion state is 
traced by the X-rays while the newly-detected 7 rays provide 
a window into particle acceleration. In early 2011, Cyg X-3 
underwent a new transition to the soft state (Kotani et al. 
2011 b), which was accompanied by transient 7 -ray emission 
detected by AGILE and Fermi (Bulgarelli et al. 201 la, b; 
Corbel et al. 2011). The transition was monitored in radio, 
7 -rays, soft and hard X-rays (§2). The Fermi / LAT detec- 
tions occurred exactly prior to and following a period of 
quenched radio emission and ultrasoft /hypersoft X-ray emis- 
sion (§3). The LAT detection accompanying a major radio 
flare casts new light on the relationship between non-thermal 
radio, X-ray and 7 -ray emission and relativistic ejection (§4). 


2 OBSERVATIONS AND DATA ANALYSIS 

Cyg X-3 is continuously monitored in 7 -ray, soft and 
hard X-ray by all-sky monitors like Fermi/ LAT, MAXI, 
RXTE / ASM and Swift/ BAT. Dedicated radio observations 
are also frequently conducted to constrain the variable ac- 
tivity of its relativistic jets. 


2.1 7-ray 

We present the results of observations of the Cygnus region 
with the Fermi Large Area Telescope (LAT; Atwood et al. 
2009) along the course of its recent 2011 active phase. The 
reduction and analysis of LAT data were performed using 
Science Tools v.9.24. LAT photons within a 15° acceptance 
cone centered on the location of Cyg X-3 were selected in 
the 100 MeV 100 GeV energy range. To minimize contam- 
ination by Earth albedo photons, 7 -ray events that have re- 
constructed directions with angles with respect to the local 
zenith > 100° have been excluded. The rocking angles were 
restricted to < 52°. Due to the potential contamination of 
the nearby pulsar PSR J2032+4127 (~ 30' from Cyg X-3), 
we used its most up to date ephemeris 1 to select the LAT 
data from its off-pulse phase intervals (removing only ~ 19 % 
of the useful exposure). The instrument response functions 
(IRFs) “P6_V11_DIFFUSE” have been used throughout this 
paper. 

Aiming to construct a light-curve on ~ daily timescales, 
we first characterized the sources within this crowded field 
of view (e.g. FLC09) by using the Fermi data over the 
two-year period with an internal source list made under 
IRFs “P6_V11_DIFFUSE” and different spectral models. A 
binned maximum likelihood spectral analysis was then per- 
formed to further constrain the spectral parameters of the 
sources within 3° (with all free parameters) and 9° (with 
only the normalization free) from Cyg X-3. The two nearby 
bright pulsars (PSR J2021+4026 and PSR J2021+3651) 
were modelled with an exponentially cut-off power-law 
model with all parameters let free to vary. We included mod- 
els for the Galactic diffuse emission (glLiem_v02) and an 
isotropic component (isotropic_iem_v 02 ) including the extra- 
galactic diffuse emission and the residual background from 
cosmic rays. Both diffuse components were renormalized for 
use under IRFs “P6.V11.DIFFUSE”. 

Once the spectral parameters of the sources in the field 
of view were fully characterized, we fixed those parameters 
and performed an unbinned maximum likelihood analysis on 
different timescales ranging from 6 hours to 4 days. In this 
procedure, the remaining free parameters are the normaliza- 
tion of the Galactic diffuse emission and the normalization of 
the power-law that is used to model Cyg X-3 (as described in 
FLC09). In time bins where Cyg X-3 is not detected, which 
we take to correspond to a test statistic, TS < 20 (4-day 
bin) or TS < 9 (shorter time bins), a 95% confidence upper 
limit is calculated. The upper limit is calculated using the 
bayesian method of Helene (1991), as implemented in the 


1 LAT pulsar timing models available here: 
https: //confluence. slac.stanford.edu/display/GLAMCOG/ 
LAT+Gamma-ray+Pulsar-|-Timing-|-Models 
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Python UpperLimits module provided with the Fermi Sci- 
enceTools, with the Cyg X-3 power law photon index fixed 
to be 2.7. 

Following the reactivation of Cyg X-3 in 7 rays 
(Corbel et al. 2011), a dedicated Target of Opportunity 
(ToO) observation was also conducted by Fermi/ LAT start- 
ing on 2011 March 24 (MJD 55644.65). However, the end of 
the flaring activity of Cyg X-3 resulted in an termination of 
the ToO observation on March 28 (MJD 55648.63). 

2.2 Radio 

Since the launch of Fermi in June 2008, we have had an 
ongoing monitoring program of Galactic binaries with the 
Owens Valley Radio Observatory (OVRO) 40 m single-dish 
telescope located in California (USA). The OVRO flux den- 
sities are measured in a single 3 GHz band centered on 15 
GHz using Dicke switching and dual-beam sky switching to 
remove atmospheric interference. The Galactic binaries sam- 
ple, including Cyg X-3, is observed with the same cadence 
and procedures used for the blazar monitoring program de- 
scribed in Richards et al. (2011). In addition, two days of 
intensive monitoring of Cyg X-3 were performed while the 
source was above the OVRO horizon in February 2011 (MJD 
55599 and 55601). An offset of 0.124 Jy has been removed 
from the OVRO flux densities to account for the effect of 
extended nearby sources (Sanchez-Sutil et al. 2008) that are 
usually resolved out by interferometers (based on a compar- 
ison with AMI data obtained simultaneously at the same 
frequency; FLC09). 

We take advantage of the large set of measurements pro- 
vided by the AMI Large array (Cambridge, UK), consisting 
of a set of eight 13 m antennas with a maximum baseline 
of ~ 120 m. AMI observations are conducted with a 6 GHz 
bandwidth centered at 15 GHz (Zwart et al. 2008). We also 
used the 11.2 GHz data of Cyg X-3 from a monitoring pro- 
gram (Trushkin et al. 2006) with the RATAN-600 telescope 
located near Zelenchurskaya vil. on the North Caucas (Rus- 
sia). Its 11.2 GHz radiometer is cooled by the cryogenic sys- 
tem up to 12 K. No frequency correction has been applied 
to the RATAN data. 


2.3 X-ray 

We used the one day average quick-look measurements in 
soft and hard X-rays from: 1) the All-Sky Monitor aboard 
the Rossi X-ray Timing Explorer (RXTE/ ASM; 1.5-12 keV; 
Levine et al. 1996), 2) the Monitor of All-sky X-ray Im- 
age (MAXI; 1.5-20 keV; Matsuoka et al. 2009) installed on 
the International Space Station, and 3) the Burst Alert 
Telescope on the Swift mission (Swift/ BAT; 15-50 keV; 
Barthelmy et al. 2005). 

We also obtained pointed RXTE observations between 
2011 February 19 (MJD 55611) and 2011 April 4 (MJD 
55655). The 46 approximately daily monitoring observations 
performed under this program (P96375) were triggered by 
Cyg X-3 entering the ultra-soft X-ray /quenched radio state, 
and we followed the evolution of the Cyg X-3 X-ray prop- 
erties as the source made its expected transition out of this 
state. We extracted the Proportional Counter Array (PCA; 
Jahoda et al. 2006) data from these observations, and the 


source evolution is characterized in terms of the hardness 
ratio and the detailed spectral properties below. The aver- 
age PCA exposure time per observation is ~ 2.7 ks. 

We produced 3-50 keV PCA energy spectra for the 46 
RXTE observations. We used the standard tools provided in 
HEASOFT v6.11 to produce source and background spectra 
and included the recommended 0.5% systematic uncertain- 
ties. We carried out the spectral fits with the XSPEC vl2 
software package. As the spectrum varies with the orbital 
modulation, we restricted our analysis to the observations 
obtained near the peak of the modulation. Table 1 lists the 
five observations for which we carried out spectral fits along 
with their orbital phases. The first of these observations oc- 
curred when the source was still in the ultra-soft state and 
the last was after the gamma-ray flare and during the decay 
of the radio flare. 

We also used the RXTE/ PCA data to characterize the 
Cyg X-3 X-ray timing properties during the five observations 
for which we reported detailed spectral properties above. 
We made light curves from the “Standard 1” data, which 
has 0.125 s time resolution and effectively no energy res- 
olution. The single energy bin includes all PCA channels, 
corresponding to an energy band of ~2-60 keV. We first pro- 
duced 0.002-4 Hz power spectra in the Leahy et al. (1983) 
normalization and then subtracted the Poisson noise level 
and divided by the total count rate to convert to the rms 
power spectra (Miyamoto et al. 1991) that are shown in Fig- 
ure 5. We fitted the power spectra with a power-law model, 
and used these fits to determine the overall fractional rms 
level for each observation. 


3 RESULTS 

3.1 Radio and X-ray emission during the 2011 
soft state 

The long-term light curves of Cyg X-3 in soft and hard X- 
rays are presented in the two top panels of Figure 1, along 
with the radio lightcurve in the bottom panel. The X-ray and 
radio states, as defined by Szostek, Zdziarski & McCollough 
(2008), are marked on top of this Figure. Up to mid-October 
2010 (~ MJD 55480), Cyg X-3 is found in its typical hard 
state. At that time, MAXI and RXTE/ ASM data show a 
rise in soft X-rays (« 2 — 5 keV) whereas Swift/BAT mea- 
sures a decline in hard (> 15 keV) X-rays, highlighting an 
initial transition to the soft/suppressed X-ray /radio state. 
The 3-5 keV count rate from the RXTE/ ASM of « 3 ctss -1 
corresponds to the pivotal value separating hard from pre- 
dominantly soft X-ray states in Cyg X-3 (see Figure 6 of 
Szostek, Zdziarski & McCollough 2008). The decrease (~ 
MJD 55500) in soft X-ray flux indicates a temporary re- 
turn to the level of the hard/soft state transition and minor 
flaring radio state. The further increase in soft X-rays above 
the pivotal ASM value starting from « MJD 55540 marks 
the full transition to the soft state, while the abrupt drop in 
hard X-ray flux on MJD 55608 (February 16, 2011) indicates 
the transition to the ultra-soft X-ray state. 

Once fully in the soft state, Cyg X-3 is again in the 
suppressed radio state (between MJD 55540 and 55570), 
and then the radio emission is variable with the pres- 
ence of flares with peak flux densities reaching ~ 0.6 Jy 
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Figure 1. Hard X-ray light curve in panel (a) from Swift /BAT (15—50 keV; green) with the soft X-ray light curves in panel (b) from 
RXTE/ASM (3—5 keV; red) and MAXI (1.5—4 keV; orange) from 2010 September 21 to 2011 June 8 . The horizontal dotted line in (a) 
and (b) highlights the emission level corresponding to the hard to soft state transition (e.g. the ASM 3-5 keV pivotal rate of 3 cts s — 1 of 
Szostek, Zdziarski & McCollough 2008). Panel (c): LAT Flux (E > 100 MeV) light curve of Cyg X-3 obtained in 4-day bins. The LAT 
fluxes (left axis) are expressed in units of 10 ' ph cm -2 s — 1 above 100 MeV. LAT upper limits are represented at the 95% confidence 
level. Panel (d): Radio lightcurve of Cyg X-3 at 15 GHz based on the OVRO, AMI and RATAN (at 11.2 GHz) data. The onset of the 
giant radio flare is indicated in the four panels by the large vertical grey band, whereas the peak of this flare is marked with the vertical 
dotted line (see Figure 3 for a zoom on this interval and Figure 2 for period I). The horizontal dashed lines in panels (a) and (d) (and 
hatched area = the 0.2 to 0.4 Jy zone) highlight the thresholds (see discussion in paper) corresponding to the detection of 7 - rav emission 
of Cyg X-3 by LAT. The labels in the top of the figure indicate the corresponding X-ray /radio state of Cyg X-3 (“Quies.”, “Sup.” and 
“Maj. fl.” stand respectively for “Quiescence”, “Suppressed” and “Major flaring” radio state). 


at 15 GHz. As opposed to the solitary major radio flares 
(> 10 Jy), which are known to follow quenched radio 
states and to decay on timescales of days, these small ra- 
dio flares (< lJy) are shorter in duration, more rapidly 
variable and may occur in clusters of groups that are 
alternating with short periods of radio emission at the 
level typical of the suppressed radio state (a good exam- 
ple can be found in Miller- Jones et al. 2009). The small 
radio flaring activity was first named by Waltman et al. 
(1994) and was not included in the classification of X-ray 
and radio states of Szostek, Zdziarski & McCollough (2008). 
Apparent superluminal expansion during two small flares, 
with peak flux up to 0.3 Jy at 15 GHz, was observed by 
Newell, Garrett & Spencer (1998). Furthermore, on several 
occasions these small flares have been observed up to the 
millimeter domain, which was also occasionally the case 


early in 2011 with inverted spectra up to 98 GHz (e.g. 
Kotani et al. 2011a; S. Trushkin private com.), suggesting 
temporary optically thick synchrotron emission. All these 
radio properties are consistent with Cyg X-3 moving in early 
2011 from the suppressed state to a state of small flaring ra- 
dio activity. 

The time span between MJD 55610 and 55640 with 
weak radio emission (down to 2 mjy) and very low hard 
X-ray emission is characteristic of the quenched radio 
state that always precedes a giant flare (Waltman et al. 
1994; Fender et al. 1997). Cyg X-3 is then in its ultra- 
soft (Szostek, Zdziarski & McCollough 2008; Koljonen et al. 
2010) X-ray spectral state characterized by strong thermal 
soft X-rays and a very weak or absent non-thermal power- 
law in hard X-rays. Optically thick radio spectra are usually 
observed in this state (Waltman et al. 1995). 
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Figure 2. Same as Figure 1 but for a time interval around period I (2010 December 30 to 2011 February 28). The soft X-rays (not shown 
for clarity) were all above the ASM pivotal value of 3 cts 1 ). The LAT fluxes in 1-day time bin are also plotted in the middle panel. 


The giant radio flare that ends the quenched state 
reached a flux of almost 20 Jy at 15 GHz (Corbel et al. 
2011). Figure 3 highlights the rising part of this flare, indi- 
cating a likely onset (defined by the gradual increase of ra- 
dio emission at the end of the quenched radio state) around 
MJD 55641.0 ± 0.5 and peak flux around MJD 55644 (also 
marked on Figure 1). The delay between onset and peak 
is consistent with observations during previous giant radio 
flares (of the order of 2 to 4 days, e.g. Waltman et al. 1995; 
Miller- Jones et al. 2004). We caution that the exact trigger 
time of this relativistic ejection event could possibly occur 
after, or even before, the radio onset we highlighted above. 
Alternatively, the detection of a non thermal X-ray compo- 
nent by RXTE / PCA on MJD 55642.0 (Fig. 4, Table 1 and 
section 4.2) and the hardening of the X-ray spectra (Fig. 3) 
may also signal the trigger time of the ejection event (in- 
dicated on Figure 3). The precise ejection date should be 
better constrained with the VLBI campaign we conducted 
during this flare (Miller- Jones et al. in prep.). 

The hard X-rays switch back on and are cor- 
related with the radio flux during major flares 
(Szostek, Zdziarski & McCollough 2008). The RXTE/PCA 
hardness ratio (Figure 3) confirms the hardening, as well 
as the appearance of a hard X-ray tail in the RXTE / PCA 
spectra (see also Figure 4 and discussion in section 4.2) 
during the onset of the flare. This hard X-ray tail was not 
present during the radio quenched period of the 2011 soft 
state (Figs. 1 and 4). The spectra during major radio flares 
typically show a non-thermal power-law in hard X-rays 
(Szostek, Zdziarski & McCollough 2008; Koljonen et al. 
2010) and optically thin radio emission (Waltman et al. 
1995). 


3.2 Flaring 7 -ray emission from Cyg X-3 

The Fermi / LAT 7 -ray light curve (0.1-100 GeV) with 4- 
day time bins is presented in panel (c) of Figure 1. Only the 
very significant data-points with test statistic (Mattox et al. 
1996) TS ^ 20 (~ 4.5 a) are plotted, otherwise, they are 
represented as upper limits. Two distinct phases of 7 -ray 
activity are found: (1) a long phase (~ 3 weeks) of 7 -ray ac- 
tivity just before the ultra-soft state (hereafter period I); ( 2 ) 
a shorter 7 -ray active phase (< 5 days) in conjunction with 
the giant radio flare (hereafter period II). No 7 -ray emis- 
sion is detected during the quenched radio state in between 
period I and II, corresponding to the ultra-soft X-ray state. 

Regarding period I (see also zoom in Fig. 2), we note 
that the 4-day bin data (from MJD 55585 to 55610) are 
consistent with a steady increase of 7 -ray emission up to a 
flux (^100 MeV) ~ 1.8 x 10 -6 ph cm -2 s -1 , just before 
transition to the quenched radio state. However, the data 
on shorter bins (Fig. 2) indicate variability on timescales as 
short as one day. 7 -ray emission in period I occurs after the 
soft X-ray emission has increased beyond the ASM pivotal 
value (~ 3 cts s -1 ), also during a significant decrease in hard 
X rays (BAT flux < ~ 0.02 cts cm ~ 2 s _1 ) and in a period of 
short and faint radio flaring, placing Cyg X-3 in the small 
flaring radio state (as discussed in section 3.1). The AGILE 
collaboration also reported detections of Cyg X-3 during pe- 
riod I (Bulgarelli et al. 201 la, b). No pointed RXTE obser- 
vations were executed in period I, and therefore we cannot 
determine if a non-thermal hard X-ray tail was present dur- 
ing that time, as we do below for period II. 

For period II, we constructed LAT light curves on 6 - 
hour, 12 -hour and 1 -day intervals to highlight the onset of 
the 7 -ray emission and to allow a comparison with the giant 
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Figure 3. Top Panel: LAT flux (E ^ 100 MeV) light curve of Cyg X-3 obtained in 1-day, 12-hour, and 6-hour bins for the period around 
the giant radio flare (2011 March 15 to 31). The orange upper limit corresponds to the 95% confidence upper limit during the Fermi 
ToO observations. Middle Panel: Radio lightcurve (same data as Figure 1, including AMI, OVRO and RATAN) of Cyg X-3 highlighting 
the onset of the giant radio flare, that gives a starting date for the radio flare of MJD ~ 55641.0 ±0.5. Bottom Panel: Evolution of the 
RXTE / PCA hardness ratio (defined as count rate in band 10-50 keV over band 3-10 keV). The error bars are smaller than the symbol 
size. The arrows indicate the time of the five RXTE spectra discussed in section 4.2 and Fig. 4. For all panels, the onset of the giant 
radio flare is indicated by the large vertical grey band (based on the increase in radio flux density) or the vertical dashed line (based on 
the appearance of a non thermal X-ray component), whereas the peak of this flare is marked with the vertical dotted line. 


radio flare. Due to the shorter integration time, only data- 
points with TS > 9 (~ 3 a) are plotted (otherwise, they 
are represented as upper limits) in Figure 3. The short 7- 
ray activity in period II occurs during the rise of the radio 
flare. The 7-ray emission turns off before the peak of radio 
emission (vertical dotted line in Figure 3). The first LAT 
detection occurs on MJD 55640.5 ± 0.5, implying a 7-ray 
turn-on consistent with the exact onset of the radio flare 
defined by the increase in radio flux density (illustrated by 
the grey area in Figure 3; see also section 3.1). Based on the 
6-hour and 12-hour bin lightcurves, the peak in 7 rays (~ 
3 to 4 x 10~ 6 ph cm -2 s -1 ) is reached on MJD 55642.75 
± 0.25, which is well before (by 1.5 day) the peak flux in 
radio. Short timescale variability might be present within 
this active phase: the LAT 1-day bin on MJD 55641.5 ± 0.5 
is consistent with a non-detection although there is no sig- 
nificant change in LAT exposure towards Cyg X-3 on this 
specific day. This would be consistent with the upper lim- 
its reported by AGILE around this time (Bulgarelli et al. 
2011c). Furthermore, we note that most of the LAT detec- 
tions occur after the onset of the non-thermal X-ray com- 
ponent (illustrated by the dashed line in Figure 3; see also 
section 3.1). Fermi conducted a dedicated ToO observation 
from MJD 55644.65 to 55648.63 with no detection of Cyg X- 
3. We derived a 95% LAT flux (^100 MeV) upper limit of 
3.1 x 10 - ' ph cm~ 2 s -1 for the ToO interval, implying a 
reduction of the 7-ray emission by more than a factor 10 on 


a timescale of a few days. Furthermore, we find no signifi- 
cant detection by the LAT at the time of the other detection 
reported by AGILE (Piano et al. 2011b) on MJD 55710. 

Integrating the LAT data separately over the two flaring 
periods results in photon indices ~ 2. 5-2. 7 consistent with 
the previously published LAT spectra (FLC09). Further- 
more, modulation of the gamma-ray emission at the orbital 
period of Cyg X-3 (FLC09) is again detected by Fermi / LAT 
in both epochs; but the short duration of the present activ- 
ity prevents a detailed analysis and will be investigated in a 
subsequent study using Pass 7 IRFs 2 which have increased 
effective area at low energies. 


4 DISCUSSION 

4.1 Conditions for detection of high energy 7 -ray 
emission from Cyg X-3 

High-energy 7-ray emission has now been reported by 
Fermi / LAT at four different epochs (Oct. /Dec. 2008, 
Jun./Jul. 2009, May 2010, March 2011). In all cases the 
LAT detections were contemporaneous, but not coincident, 
with the hard X-ray lightcurve crossing a threshold level, 

2 Information available here: http://fermi.gsfc.nasa.gov/ssc/ 

data/analysis/documentation/Pass7_usage.html 
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Spectral and timing parameters for the 5 RXTE observations of Cyg X-3 


Parameter 

Spectrum 1 

Spectrum 2 

Spectrum 3 

Spectrum 4 

Spectrum 5 

ObsID 

29-00 

32-00 

33-00 

36-00 

40-00 

MJD-55000 

639.0 

642.0 

643.0 

646.0 

649.8 

Orbital Phase 

0.46-0.65 

0.51-0.64 

0.41-0.60 

0.46-0.65 

0.43-0.60 

Spectral 

W*s 

0 1 5+ 0 - 05 
U'-LO-O.Ol 

0.28ig;“ 

°-018±°;“o5 

0.061 j^Qog 

0.034 ±0.016 

kT bb (eV) 

319±9i 4 

291+ 1 ® 9 

247+ 12 

2851 1 ® 

257+J® 

lnt/lh 

<9 x 10 “ 4 

0 01 3+ 0 - 002 
U.UTO_0 003 

0 22 +0 01 
u -^ -0.03 

0 44 + 0 - 26 
0.12 

0 996 -1 " 0 004 

u.y»u _ 0 500 

t p 

o g+io 

°-°-o.i 

7 7 + 1 - 7 
'• ' - 1.0 

0 026 +0 008 
U.UZD _ 0 ooi 

0 049”*” 0 009 
u.u^y _ 0 019 

0 Q3Q + 0 069 
U.UoU_Q 023 

Tinj 

4 5 +0 ' 5 

4 o + 00 

5 0 +0 -° 
°- U -0.5 

tr 0 + O .2 
0.3 

c O+0.2 
^ -0.2 

x 2 + 

78/62 

57/62 

52/63 

61/62 

46/62 

Timing 

a 

1.97 ±0.20 

o 9+2.1 

° -0.9 

2.7 ±0.4 

2 29" 1 ” 0,13 

Z ‘ Zy -0.14 

2 14+ 0 - 24 
-0.27 

r.m.s. (%) 

o 91 + 0.02 

Z ’ Z1 -0.14 

0 7 C 5 + 0 - 01 
u - ' °-0.23 

1 yq + 0.02 

± - '°-0.18 

o >70 + 0.13 
0.29 

2 09”*” 012 

z.u»_o 30 


Table 1. Main spectral parameters (according to their definitions in eqpair; see 
Hjalmarsdotter et al. 2009) from the best fitting models of the five RXTE spectra discussed 
in section 4 and Fig. 4. The eqpair parameters: l b /l s defines the ratio of the hard to soft com- 
pactness, kT bb is the temperature of the inner edge of the accretion disk, Z n t/£h corresponds to 
the fraction of power supplied to the energetic particles which goes into accelerating non-thermal 
particles, r p is the thermal plasma optical depth and T7,j is the index of the power-law of the 
accelerated non thermal electrons. For the timing section in the bottom of the table, we report 
the power-law index a used in the fitting of the power density spectra (see Figure 5) along with 
the rms amplitude (in %) of the variability in the 0.002-0.1 Hz range. The quoted uncertainties 
correspond to 90% confidence. 


i.e. Swift / BAT count rates of 0.02 cts cm -2 s ” 1 (if only for 
a day in 2010, Williams et al. 2011). In period I, the cross- 
ing proceeded from higher to lower flux, whereas in period 
II it was in the opposite direction. The RXTE / ASM flux in 
the 3-5 keV band was also always at or above the pivotal 
value of 3ctss _1 that separates the canonical hard and soft 
X-ray states of Cyg X-3 (Szostek, Zdziarski & McCollough 
2008; Koljonen et al. 2010). These conditions on the X-ray 
fluxes appear necessary to detect 7 -ray emission (see also 
Piano et al. 2011a for the AGILE detections). 

When its X-ray spectrum is predominantly soft, Cyg X- 
3 can be found in a suppressed, quenched, small or major 
flaring state depending on the behavior in radio. The 2008-9 
detections and the new LAT detection in period II occurred 
during major flaring activity, with radio fluxes in excess of 
a few Jy and optically thin spectra. The new LAT detection 
in period I occurred at a time when Cyg X-3 was most likely 
moving through a series of small radio flares, with 15 GHz 
flux densities « 0. 1-0.6 Jy. The radio and X-ray conditions 
during the brief 7 -ray detection of May 2010 (Williams et al. 
2011) appear similar to those of period I. We do not detect 
7 -ray emission during the quenched radio state, when the 
radio emission is very low and the X-ray spectrum is ul- 
trasoft. However, we can not exclude that very faint 7 -ray 
emission (at least much below the sensitivity of LAT) could 
eventually be present in this state, as the radio lightcurve 
(Fig. 1) indicate some variability from the jets. We also note 
that 7 rays are not present in the mid to later stages of pe- 
riod II, even though the radio emission is still very bright 
(> 1 Jy). Slowly decaying jets are not associated with 7 -ray 
emission, and consequently a rapid rise in radio emission 


(possibly caused by strong shocks, see section 4.3), up to at 
least ~ 0.2 to 0.4 Jy, may therefore be required. 

Hence, the level of 7 -ray emission depends also on the 
presence of radio emission at a level greater than ~ 0.2 to 
0.4 Jy, reinforcing the association of 7 -ray activity with the 
significant emission from the relativistic jets. In other words, 
particle acceleration to very high energies does not happen 
exclusively during rare major radio flares but a relatively low 
level radio activity suffices. To summarize, three conditions 
seem to be required in order to detect significant high energy 
7 -ray emission from Cyg X-3: (1) a high level of soft X-ray 
emission (the 3-5 keV RXTE / ASM value above 3ctss _1 , i.e. 
Cyg X-3 needs to be in the soft state), (2) a low level of hard 
X-ray emission (Swift / BAT below 0.02 cts cm ' 2 s'” 1 ), and 
(3) the presence of significant emission with rapid variation 
from active relativistic jets (with 15 GHz radio flux above 
« 0.2 to 0.4 Jy). However, 7 -ray emission can temporarily 
be seen at lower radio fluxes (e.g. onset of emission during 
epoch II), and therefore a comprehensive picture taking into 
account possible delays or binning will be examined in future 
work. Furthermore, it is unclear at present how the shape 
of the radio spectrum (optically thin or thick synchrotron 
emission) relates to the 7 -ray emission. 

4.2 A non thermal X-ray component during the 
giant radio flare 

Giant radio flares start in the ultrasoft X-ray state, which 
also shows a rising non-thermal component beyond 15 
keV while the soft thermal emission remains steady (e.g. 
Szostek, Zdziarski & McCollough 2008). The RXTE / PCA 
spectra during the transition (Figure 4) confirm the appear- 
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Energy (keV) Frequency (Hz) 


Figure 4. RXTE / PCA spectra during the soft to hard state 
transition associated with the giant radio flare of Cyg X-3 in 2011. 
They correspond to the periods: before the onset of the radio and 
7 -ray flare (spectrum 1 on MJD 55639), (2) the onset of the 7 -ray 
flare (spectrum 2 on MJD 55642), (3) the peak of the 7 -ray flare 
(spectrum 3 on MJD 55643), and (4) after the peak of the radio 
flare (spectrum 4 and 5). For illustrative purposes, adjacent bins 
have been grouped until they reach a significant detection of at 
least 3cr, but with a maximum of 10 bins combined. 

ance of a hard X-ray tail during the onset of the flare in 
period II. The March 2011 giant radio flare shows that the 
simultaneous rise in non-thermal radio and hard X-ray com- 
ponents also involves high-energy 7-ray emission. One of the 
RXTE / PCA observations occurred on MJD 55643 (spec- 
trum 3 in Figure 4 and Table 1) while the source was being 
detected in the 7-ray band by LAT. Due to the complicated 
shape of the X-ray spectrum, the extrapolation from the 
X-ray to the 7-ray band is not straightforward. However, 
a rough estimate can be made by re-fitting the 10-50 keV 
portion of the PCA spectrum with a power-law. The value 
of the best fit photon index is F = 3.4, and an extrapolation 
to >100 MeV predicts a flux of 1.0 x ICC 11 ph cm -2 s^ 1 . 
This is several orders of magnitude lower than the measured 
LAT flux, indicating that the X rays and 7 rays cannot be 
part of the same power-law (even if the two emission regimes 
could still be connected) and that a spectral break is located 
between ~ 100 keV and 100 MeV (see also Zdziarski et al. 
2011). If the hard X rays and 7 rays are due to the same ra- 
diative process and generated at the same location then the 
orbital modulation at both frequencies should be in phase. 

It is well known that Cyg X-3 has a highly complex 
X-ray spectrum, and, as expected, simple models such as 
a power-law, a power-law with a cutoff, or thermal Comp- 
tonization (all with simple absorption) provide very poor 
fits. Thus, we used a model that has previously been used for 
Cyg X-3 (Szostek & Zdziarski 2008; Hjalmarsdotter et al. 
2009), which is based on the eqpair hybrid thermal/non- 
thermal Comptonization model (Coppi 1992, 1999). In ad- 
dition to eqpair, we included absorption with partial cov- 
ering with two separate values of Ah and iron features (a 
broad line, a narrow line, and an edge) as described in 
Hjalmarsdotter et al. (2009). As shown by the \ 2 values in 
Table 1, this model provided acceptable fits. 


Figure 5. RXTE / PCA power density spectra during the soft to 
hard state transition (same dates and colors as in Figure 4). 


The eqpair parameters from the fits to the X-ray 
spectra are given in Table 1. Whereas the purpose of this 
work is not to do detailed spectral modeling (e.g. see 
Szostek & Zdziarski 2008; Hjalmarsdotter et al. 2009), this 
information can be used to provide some insights into the 
evolution of the X-ray spectra during the major radio flare 
(period II). However, as we do not have X-ray data below 
3 keV, several degeneracies within the eqpair model (e.g. 
the hydrogen column density and the inner accretion disk 
temperature) cannot be removed. Whereas the model pro- 
vides good fits to the data, the values reported in Table 

1 should not be over-interpreted. The solutions we obtained 
are in the regime of low seed photon temperature fcTbb, very 
soft electron injection spectrum Fi n j and low plasma optical 
depth r p . Within this regime, one trend that is notable is 
the gradual increase in the l n t/lh parameter. This appears 
to be an indication of the increasing importance of non- 
thermal Comptonization in the spectrum as Cyg X-3 leaves 
the ultra-soft state. It is not clear how this is related to the 
gamma-ray emission since the LAT detected Cyg X-3 only 
during the acquisition of spectrum 3, and this is not the 
spectrum with the largest non-thermal X-ray contribution 
(unless delayed emission is involved). 

The power density spectra (PDS) corresponding to the 
five RXTE / PCA observations discussed above are plotted 
in Fig. 5. The PDS at frequencies higher than 10~ 3 Hz 
are well described by a power-law with an index of ~ - 

2 (see the timing parameters in Table 1). No signal is de- 
tected above 0.1 Hz. The slope and rms noise levels of the 
PDS are similar to previous studies (Choudhury et al. 2004; 
Axelsson, Larsson & Hjalmarsdotter 2009). One noticeable 
difference can be seen in the power spectrum # 2 with a 
lower noise level. This observation corresponds to the onset 
of non-thermal Comptonization as discussed in the previous 
paragraph. This lower variability may possibly be associated 
with a quiet accretion disk between the onset of the radio 
flare and the formation of the corona producing the hard 
X-ray emission. 
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4.3 Connecting the relativistic jets and the high 
energy emission in Cyg X-3 

A coherent picture of the link between accretion, ejection 
and the non-thermal emission we observe has yet to emerge. 
High levels of soft X-rays together with significant radio 
emission are related to the 7 -ray activity but the relative 
timings at the various wavelengths remain confusing: par- 
ticle acceleration and cooling may lead to differing lags 
and peaks; correlations may also be blurred when multi- 
ple flares superpose. Rising soft X-ray emission is thought 
to accompany an increased accretion rate (perhaps due to 
variations in mass loss rate from the Wolf-Rayet companion 
star, e.g. Gies et al. 2003) and a decreasing inner disk ra- 
dius, eventually quenching the jet in the ultrasoft state (e.g. 
Hjalmarsdotter et al. 2009). 

VLBI observations have shown that the variations in 
radio flux during major flares come from the resolved jet 
on milliarcsecond (7 a.u.) scales and not from the core 
(Tudose et al. 2010), raising the intriguing possibility that 
the non-thermal emission region is outside the binary sys- 
tem during major flares. However, there is no evidence yet 
that the radio emission is detached from the core during 
the rising phase of a major flare (when gamma-ray emis- 
sion was detected). Locating the high energy electrons very 
far away (e.g. via VLBI observations) would place stringent 
constraints on jet parameters if the observed 7 -ray modu- 
lation is due to inverse Compton upscattering of photons 
from the Wolf-Rayet star (Dubus, Cerutti & Henri 2010; 
Sitarek & Bednarek 2011; Zdziarski et al. 2011). 

A possible scenario is that the non-thermal emis- 
sion is related to shocks forming at various distances 
along the jet, as previously suggested by modeling of 
the radio activity (Lindfors et al. 2007; Miller- Jones et al. 
2009). Transitions in/out of the ultrasoft X-ray state 
then signal a decrease/increase in jet efficiency with the 
non-thermal region moving in/out. Gamma-ray emission 
may be most efficient at some “sweet-spot” distance 
bounded by strong pair production on thermal X-rays 
(Cerutti et al. 2011; Sitarek & Bednarek 2011) and a de- 
clining seed photon density for inverse Compton scattering 
(Dubus, Cerutti & Henri 2010; Zdziarski et al. 2011). De- 
tections prior to and after the quenched state would be due 
to the shock moving through this region as the jet turns 
off/on. 

The application of the shock-in-jet model to Cyg X-3 
strongly suggests that the faster, weaker radio flares (like 
in period I) occur closer to the core, whereas the brighter 
radio flares (like in period II) occur further downstream 
(Miller- Jones et al. 2009). A shock closer to the core dur- 
ing period I than during period II is also consistent with the 
brighter 7 -ray emission that is observed in period I, assum- 
ing a jet with constant speed. The energy density in seed 
photons decreases if shocks occur far downstream, reduc- 
ing the inverse Compton luminosity in period II. An addi- 
tional signature should be a stronger 7 -ray modulation when 
shocks occur close to the core, hence close to the Wolf-Rayet 
companion star. 
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